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Abstract 

In light of the recent work by Sahni et al, Harbola, and Gaudom and Burke, the question of 
mapping from an excited-state density of a many-electron interacting system to the potential of 
the related non-interacting system is analyzed. To do so, we investigate the Levy-Nagy criterion 
quantitatively for several excited-states. Our work indicates that Levy-Nagy criterion may fix the 
density to potential map uniquely. 
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The question of whether there exists a mapping from an excited-state density p(r) to 
a potential v{r) is central to performing density-functional calculations for excited states. 



le issue has been addressed recently in a series of papers by Sahni et al. Ill, Harbola 

n n n 

I, and Gaudoin and Burke [3|. In the work of [1] and [2], it was shown that a given 
ground- or excited-state density can be generated as a noninteracting-system density by a 



configuration o, 
virial theorem 



one's choice. Sahni et al. obtained the potentials using the differential 

sj. It is 



^, whereas Harbola did so using the constrained-search approach 
clear that because of the multiplicity of potentials that could lead to a given density, one 
needs an additional condition for mapping the density p(r) to a unique potential v{r). For 



the ground-state density the Hohenberg-Kohn (HK) theorem p fixes the Kohn-Sham (KS) 



system [7| uniquely - it is that system where the lowest energy orbitals are occupied. For 
the excited-state density a different criterion is needed. The issue of how to choose a unique 
potential (KS system) for a given excited-state density has been addressed earlier by Levy 
and Nagy (LN) 8]. They have proposed a qualitative criterion, discussed below, for doing 
so. However, it has not been investigated quantitatively. 



Furthering the work of [![ and 



, Gaudoin and Burke [3| have shown that even with a 
fixed configuration, one can reproduce an excited-state density from more than one poten- 
tials. They have worked with non-interacting fermions and have generated these potentials 
using the inverse linear response of a system and have related the multiplicity of these po- 
tentials to a property of the linear response kernel. Based on this they have made two 
observations: (i) that one has to go beyond the HK 3| theorem to understand the multiplic- 
ity of potentials with the same density, and (ii) that the mapping p(r) v{r) is not unique. 
However, as pointed out above, the uniqueness between a density and a potential exists 
only for the ground-states. Thus the existence of more than one potential for excited-states 
densities is not excluded by the HK theorem. And precisely for this reason, an additional 
condition is needed to identify one particular system as the KS system representing an 
excited-state. 



The interesting and thought provoking results of references 



m 



have prompted us to 
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pursue the matter of density-to-potential mapping for excited-states further. Our investiga- 
tions in this direction form the contents of this paper. We show: (i) that the constrained- 
search approach ^] itself is capable of generating all potentials for noninteracting systems 
ivingthe same density, thereby establishing once again the importance of this approach 
1 ^ 3 excited-state density-functional formalism; (ii) in the examples taken, when 
the criterion proposed by LN jsl is applied to different noninteracting systems corresponding 



to a given density, it correctly identifies the system that should represent the excited- state 
density; and (iii) the conclusion of Gaudoin and Burke about the "lack of HK theorem 
for excited-states" has been arrived at by taking into consideration only the excited-state 
density p(r) and therefore do not apply to the Levy-Nagy formalism. In the following, after 
describing briefly the LN formulation of density-functional theory (DFT) for excited states, 
we give examples of how its application gives a unique p(r) — >■ v{r) map. 

Like for the ground-state theory, the LN formulation provides a variational DFT 
approach for the kth excited state of an A^-electron interacting system by defining a unique 
universal functional Ffc[p, po] such that the energy and the density Pfc(r) of this state are 
given by 

Ek = min | j d^r t;cxt(r) p(r) + Fk[p, po] | 

= j dV Ve^tir) pfc(r) + Fk[pk,Po] (1) 
Here fext(r) is the external potential, po(r) is the ground-state density of this system. Due 



to the HK 
functional 



theorem, fext is a unique functional of po- In the definition of the bi-density 



F,[p,po] = rniri (vl/|f + Kel^) , (2) 

{{•^\^j)=0, j<k} 



the A^-electron trial wave function \1/ belongs to the space which is orthogonal to the space 
spanned by all lower jth-state functions of the system, j < k. T and Vee are the A^-electron 
operators of the kinetic and electron-electron interaction energies. In this formulation, the 
mapping from the kth excited-state density Pfc(r) to the corresponding wave function of 
the system follows from Eq. ^ after inserting there p = Pk — the minimizer in Eq. (Q), 



because then is the minimizer in Eq. Since the lower states (j < k) in Eq. (j21) 
are determined by fext[Po]) a functional of po, it is clear that the ground-state density po 
plays an important role in this DFT for excited states. 

For each original, interacting system, one can introduce a corresponding noninteracting 
system such that their densities are the same. In the constrained-search approach, this is 
done by minimizing the expectation value (•I'lTl^), where $ is now a single Slater determi- 
nant of one-electron spin orbitals that gives the density of interest. However, many different 
noninteracting systems (potentials) can be related with the given excited-state density pfc(i") 
of the interacting system. Of the many $'s (many systems) that may give the same density 
Pk, a unique one, the KS system, is chosen by comparing the ground-state densities of the 
corresponding noninteracting systems and the true ground-state density po, and checking if 
the LN criterion for identifying the Kohn-Sham system for an excited-state is satisfied. Thus, 
let in a particular noninteracting system [characterized by its potential v{r)] the density of its 
mth state, pj^(r), be the same as Pfe(r). Its ground-state density will be denoted accordingly 
as Po(r)- Then the KS system connected with p^ is identified among the above noninteract- 
ing svstems as the one whose Po(r) resembles po(i') most closely in a least-squares sense. The 
LN 8[ criterion intuitively defines the KS system consistent with the adiabatic connection 
to the kth excited-state of the interacting system (characterized in DFT by Fk[p, po])- What 
it means is if the electron-electron interaction in an interacting system is turned off slowly, 
keeping the excited-state density unchanged, the corresponding ground-state density of the 
resulting system will remain close to the true ground-state density of the interacting sys- 
tem. Thus of the many noninteracting systems that give the same excited-state density, the 
one whose ground-state density remains closest to the true ground-state density of a given 
system is identified as the KS system representing the excited-state of that system. The 
noninteracting system so chosen should best resemble the true system because, within the 
constraint of the equality of their excited-state density, their ground-state densities match 
most closely. This should also make their external potential resemble each other by the HK 
theorem. We reiterate that in general a particular density can be generated by a multitude 



4 



of potentials; and for each potential the associated energy functional is different |2| due 
to the difference in the noninteracting kinetic energy of each system. However, to keep the 
structure of these functionals and the corresponding potentials simple, it is important that 
we have a criterion to choose one particular system. In this connection we note that for 
the ground-state densities of noninteracting electrons too, there exist 0] more than one 
noninteracting systems that give the same density po- However, the one where the lowest 
energy orbitals are occupied, i.e. Po(^) — Polr), is the chosen KS system, and it is unique 
due to the HK theorem. For the excited-states the uniqueness should be provided by the 
LN criterion. But, before applying it in practice, this qualitative LN [8J criterion that "Po(r) 
resembles po(r) most closely in a least-squares sense" needs to be transformed into some 
quantitative form. 

One of the ways that the difference between two densities pa and pf, can be characterized 
quantitatively is by the squared distance in the functional space 

A[p„,p,]= [ dV(p,(r)-p,(r))'. (3) 

We propose to consider the value of A[pq, po] as representing the least-squares deviation of the 
density Pq from the density po. Then, applying the LN criterion, the noninteracting system 
having the smallest A would be chosen as the KS system. We are going to demonstrate on 
examples that the proposed quantitative version of the LN criterion chooses the KS system 
in agreement with intuitive expectations in the cases considered. What it means is if the 
excitation corresponds (i) to a fixed external potential with several configurations(Sahni et 



al. m and Harbola 2]) or (ii) to a particular configuration with different potentials (Gaudoin 
and Burke ^), then the minimum deviation occurs only for the true configuration / for the 
exact potential. It is also shown that if some particular, different quantitative measure of the 
distance between densities is used in the LN criterion, it may lead to erroneous choice of the 
KS system. Thus definition of an adequate quantitative form of the criterion is important 
and needs verification. 

We now discuss the case of the excited-state density of a noninteracting Fermionic systems 
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where the LN criterion is very transparent. This is because in this case the difference 
between the ground-state density of the KS system representing an excited-state and the 
true ground-state density will be zero. Thus it is easily shown that the LN criterion fixes 
the KS system to be the true system: Suppose an excited-state density is produced by 
two different potentials jsj. Each of these potentials (systems) has a unique ground-state 
by the HK theorem. Thus only that particular system which has the same ground-state 
density as the system under consideration can truly represent the original system. In the 
interacting system, the condition of the ground-state densities being the same is replaced 
by the condition of the proximity of the densities in the least-square sense. We note that 
in the work of Gaudoin and Burke ^ who have analyzed non-interacting systems only, 
this comparison with the ground-state density is not made. Hence they raise a question 
which potential should be chosen from the many available. It is clear from the discussions 
above that this should be determined by the comparison of the ground-state densities of the 
excited-states systems (which are connected with the given p^) with the true ground-state 
density of the system under consideration. In the following we show through examples that 
this gives a KS noninter acting system which is consistent with the original system not only 
in terms of their ground-states, but also in terms of the configuration of the excited-states. 

To construct, for an assumed electronic configuration, the potentials leading to a given 
de„s.y thjo.|co„st.a.ed-sea.cK approach Q. we have employed the procedu. of Zhao 
and Parr Il2| . It produces the potential f (r) of a noninteracting system by making 
the value of the system kinetic energy Xli^i ~ IV^lV'j) stationary with the constraint 

that the system density ^j|'?/'j(r) P equals the given density Pkij). Here, V'i(r) is the 
space orbital, — its occupation number, Xli^i "^i ~ configuration connected with a 

single Slater determinant can be represented by a proper choice of Ui from the numbers 0,1,2. 
Depending on the starting potential used to initialize this procedure, various potentials can 
be generated from the input pk for each assumed configuration. Thus this procedure is 
capable of generating all possible systems (potentials) that reproduce the density on hand. 
Evidently, the determined different potentials give different ground-state densities (by the 
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Hohenberg-Kohn theorem [6|); the one which is the closest to the true ground-state density 
singles out the KS potential of the excited-state DFT. 
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FIG. 1: Two potentials (lower panel) yielding the same excited state density (upper panel) for 
ls^2s^ state of a model Li atom. Note that the x-axis scale in the upper and the lower panel is 
different. 
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FIG. 2: Shown are the ground-state densities for the two potentials vi and V2 of Fig. ^ 



As the first example we take a model Li atom: = 3 noninteracting electrons moving 
in the potential v{r) = —Z/r, Z = 3, resulting in the hydrogen-like orbitals. Atomic units 
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are used throughout. The k = ls^2s^ excited-state configuration of the model system is 
considered. The potential fi(r) = v{r) and the alternative potential V2{r) along with the 
excited-state density p(r) = Pfc(r) are shown in Fig. ^ Both potentials are generated using 



the Zhao-Parr 
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Il2i | method with the excited-state density pk as the input. Since the 
potentials Vi and V2 are different, their ground-state densities can not be the same (by the 
HK theorem jQ]). Whereas the ground-state density corresponding to f 1 = —Z/r is the true 
ground-state density of the system, that corresponding to f 2 should be different — they are 
shown in Fig. El Indeed, the two densities are dissimilar. If f 2 also were to represent the 
KS system connected with the same excited-state density, the difference in the ground-state 
densities for these two potentials should be zero. 




r (a.u.) 

FIG. 3: Two potentials (lower panel) yielding the same excited state density (upper panel) for 
the ls^2s^ ^/S state of the Li atom. Note that the x-axis scale in the upper and the lower panel is 
different. 

As the second example we take the true Li atom: = 3 interacting electrons moving in 
the external potential fext('") = ~Z/r, Z = 3, and consider its k = ls^2s^^S' excited-state 
configuration. The density Pfc(r) of this state is represented by the density calculated self- 
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FIG. 4: Shown are the noninteracting ground-state densities for the two potentials vi and V2 of 
Fig. 121 along with the interacting (exact) one. 

consistently using the exchange-only Harbola-Sahni(HS) 2| potential for this configuration. 
This density is a good approximation to the exact solution and very close [ij, Q to the 
Hartree-Fock density of this state. Shown in Fig. El are the two potentials fi(r) and V2{r) 
reproducing the same density pfc(r) as the ls^2s^ excited-state densities of noninteracting 
systems; vi coincides with the HS effective potential, fi(r) = —Z/r + vu{r) + t>^^(r). The 
constrained-search procedure, mentioned earlier, was employed to generate both Vi{r) and 
f2(r). Although the excited-state densities of the two potentials are the same, the ground- 
state densities of these potentials are different — that is what discriminates between the two 
potentials. In Fig. 01 we plot the noninteracting ground-state densities Po(^) ^'^^ potentials 



vi and V = V2, along with the interacting ground-state density po{r). T 



le latter is 



obtained in the same approximation as applied for pk, namely with the HS 1^ potential, 
now for the ground-state ls^2s configuration. As expected, the three densities are different. 
However, the noninteracting ground-state density produced by Vi is quite similar to the 
"exact" HS ground-state density of Li. On the other hand, that corresponding to V2 is very 
different from the "exact" one. We introduce also the potential vs{r) (not shown), which 
reproduces pfc(r) as the ground-state density of a noninteracting system. This is unique 
according to the HK theorem. Thus the density shown in the upper panel of Fig. El can also 
be labeled with in analogy with Fig. |31 To apply the LN criterion, the squared distance, 
Eq. 0, between ground-state densities is evaluated, giving A[pQ,po] = 0.111,1.467,0.460 
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for V = vi,V2,V3, respectively. Thus, according to this criterion, the KS system connected 
with the kth excited state of the Li atom is given by the potential Vi. This result confirms 
our intuitive expectation. 




-V, (X) 



-V,(x) 




FIG. 5: Two potentials (middle panel) yielding the same excited-state density (upper panel) along 
with their corresponding ground-state densities (lower panel) for an excited state of the three- 
electron ID infinitely deep well model system. 



As the third example we take the He atom and consider its k = ls2s excited-state, with 

n n 

the density pk taken from [16||. This example was examined, in fact, in the previous work [2| 
of one of us. Two potentials were obtained that reproduce the density pk- f i(r) as the density 
of the ls2s excited state of a noninteracting system, and f 2(r) as the density of the Is^ ground 
state. To make use of the LN criterion, the quantity A[pq, pq] = dr {pQ^r) — Po{r)Y (for 
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spherical densities) was defined, and evaluated with p„(r). the ground-state density of the 

true He atom, taken from [17]. On the basis of the resuhs A = 0.273,0.140 for v = vi,V2, 
respect... the author conceded that the L.er.ter.nB..ht not. proper for finding 
the KS potential of the excited-state DFT. However, if the distance between the ground- 
state densities is evaluated applying the definition (proposed in the present paper), the 
result A = 0.086,0.091 for v = fi,f2, respectively, is obtained. Thus on the basis of the 
LN criterion, we find that indeed the ls2s ^S* state of He is properly represented by the KS 
system with the potential vi that reproduces pk as the density of its ls2s configuration, in 
agreement with our intuition. 

We note that although the considered examples seem to suggest that comparing A as 
given by Eq. 3 gives KS system in accordance with the excited-states, Eq. 3 is not the only 
way of quantifying the LN criterion. Better and more discriminating criteria may exist and 
should be looked for, particularly because of the limited number of examples considered in 
this paper and very small difference that is there for the He atom. But the emphasis in this 
paper is on showing how an additional condition of comparing the ground-state densities 
may lead to a proper choice of the KS system and that is shown amply by our work. Further 
work along these lines is in progress and will be reported in the future. 

The above arguments apply equally well to the one- dimensional (ID) case considered 
by Gaudoin and Burke P|. Using another conventional approach — the van-Leeuwen and 



Baerends 18| method — for obtaining the noninteracting-system potentials, we have repro- 
duced not only the results of Gaudoin and Burke but have also done many other calculations. 
The example we give in this paper is for the following ID model system: = 3 noninter- 
acting electrons in an infinitely deep box of unit length, in the excited state obtained by 
putting one electron in the lowest-energy state and two electrons in the second-lowest one. 
The excited-state density (pe) and the corresponding two potentials that reproduce this den- 
sity in the same configuration are shown in Fig. El Also shown are the ground state densities 
Po^ and po^ corresponding to these potentials. Again the ground-state density given by the 
potential V2{x) is not the same as that given by vi{x). Thus V2{x) can not represent the 
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excited-state density of the model. 

To conclude, we have shown that if the ground-state density is known, then the Levy-Nagy 
criterion may provide a proper map from an excited-state density to the KS potential, p(r) — »• 
f (r), provided the closeness of two densities is adequately quantified, e.g., as in Eq. Q. The 
criterion is exact for systems of non-interacting fermions. For interacting electron systems, 
on the other hand, there could be other ways of defining this closeness quantitatively, but 
their effect should be verified (an example of failed definition was discussed). Our focus 
here is not on various definitions but rather how the proposed measure leads to a map from 
an excited-state density to the corresponding Kohn-Sham potential that is consistent with 
the configuration of excitation in the known cases. Thus the excited-state energy can be 
expressed in terms of the density corresponding to this state and Kohn-Sham calculation 
can be done following the Levy-Nagy formulation. 
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